An Apparatus and Method For Providing a 
Flow Measurement Compensated for Entrained Gas 
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Technical Field 

I 5 This invention relates to an apparatus for measuring a flow having entrained gas 

therein, and more particularly to an apparatus that measures the speed of sound propagating 
through the flow to determine the gas volume fraction of the gas in the process flow and 
compensating the output measurement of a flow meter (e.g., a volumetric flow meter and a 
consistency meter) for entrained gas. 

20 

Background Art 

The present invention provides an apparatus and method of measuring volumetric 
flow rate and gas volume fraction in slurries used in the paper and pulp industries and in 
other industries. Slurries commonly used in the paper and pulp industry are mostly water 

25 and typically contain between 1% and 10% pulp content by mass. Monitoring the flow rate 
and consistency of the slurry can lead to improved quality and efficiency of the paper 
production process. 

Processes run in the paper and pulp industry can often, either intentionally or 
unintentionally, entrain gas/air. Typically, this entrained air results in measurement errors 

30 in process monitoring equipment such as volumetric flow measurements and consistency 
meters. 
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Industry estimates indicate that entrained air levels of 2-4% are common. Since 
most process flow monitors are unable to distinguish between air and liquid, interpreting 
their output as liquid flow rates would result in a overestimate of the liquid by the 
volumetric flow rate of the air present at the measurement location. Similarly, for the void 
fraction of the air within the pipe can cause errors in consistency measurements. 

Thus, providing a method and apparatus for measuring entrained air in paper and 
pulp slurries would provide several benefits. Firstly, it would provide a means to screen the 
output of process instrumentation. Secondly, in addition to screening the measurements, an 
accurate measurement of the entrained air would provide a means to correct the output of 
volumetric flow meters and consistency meters. Thirdly, monitoring variations in the 
amount of entrained air in a given process could be indicative of process anomalies, such a 
worn bushing or cavitating pumps and/or valves. 

Multiphase process flow rate is a critical process control parameter for the paper and 
pulp industry. Knowing the amounts of liquid, solids and entrained gases flowing in 
process lines is key to optimizing the overall the papermaking process (Matula, 2000). 
Unfortunately, significant challenges remain in the achieving accurate, reliable, and 
economical monitoring of multiphase flow rates of paper and pulp slurries. Reliability 
challenges arise due the corrosive and erosive properties of the slurry. Accuracy challenges 
stem from the multiphase nature of the slurries. Economical challenges arise from the need 
to reduce total life time cost of flow measurement, considering installation and maintenance 
costs in addition to the initial cost of the equipment. 

Currently, there is an unmet need for multiphase flow measurement in the paper and 
pulp industry. Real time flow measurement is typical restricted to monitoring the total 
volumetric flow rate in a process line without providing information on the composition of 
the process mixture. For example, electromagnetic flow meters are the most widely used 
flow meters in the paper and pulp industry, however they provide no indication of presence 
of entrained air, with its presence resulting in an over prediction of the volumetric flow of 
process fluid by the amount of air entrained. Consistency meter provide a measurement of 
the percentage of solids within the process, however this technology remains more of an art 
than a science. Furthermore, although entrained air is known to have a large, often 
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deleterious, impact on the paper making process, instrumentation is currently not available 
to provide this measurement on a real time basis. 

The present invention an accurate, reliable multiphase flow measurement in the 
paper and pulp industry. 

5 In one embodiment of the present invention, the apparatus and method improves the 

determination of consistency of paper and pulp slurries. Consistency refers to the mass 
fraction of pulp contained in water and pulp slurries used in the paper making process. 
Consistency measurements are critical in the optimization of the paper making process. 
Currently, many companies produce consistency meters employing various technology to 
10 serve the paper and pulp industry. Unfortunately, accurate and reliable measurement of 

consistency remains an elusive objective. Typically, interpreting the output of a consistency 
meter in terms of actual consistency is more of an art than a science. 

Of the various types of consistency meters on the market, microwave based meters 
may represent the best the solution for many applications. One such microwave-based 
15 consistency meter is manufactured by Toshiba. Microwave consistency meters essentially 
measure speed or velocity the microwave signal propagates through the medium being 
measured. For example, the speed of the microwave signal through water is approximately 
0.1 time the speed of light in a vacuum (c), through air is approximately 1.0 times the speed 
of light in a vacuum, and through fiber (or pulp) is approximately 0.6 times the speed of 
20 light in a vacuum. 

The velocity of the microwave signal propagating through the paper pulp slurry is 
measure by the conductive effects of the slurry, in accordance with the following equation: 
V-c*sqrt(E) 

Where V is the velocity of the microwave signal propagating through the slurry, c is 
25 the speed of light in a vacuum, and E is the relative conductivity of the material. Typical 
values of relative conductivity for material comprising a paper/pulp slurry, for example, are: 
Water relative conductivity = 80; 
Air relative conductivity = 1; and 
Fiber relative conductivity = 3. 
30 These meters typically work well in the absence of entrained air. With entrained air 

present, the air displaces water and looks like additional pulp fiber to the microwave meter. 
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Thus, uncertainity in the amount of entrained air translates directly into uncertainty in 
consistency. 

Summary of the Invention 

Objects of the present invention include an apparatus having sensor for determining 
the speed of sound propagating within a pipe for determining the gas volume fraction of a 
process flow to correct the output of a meter for entrained gas, such as a volumetric flow 
meter and a consistency meter. 

According to the present invention, an apparatus for measuring a parameter of a 
process flow flowing within a pipe includes a first meter portion and a second meter 
portion. The first meter portion provides a meter measurement signal indicative of a 
parameter of the flow propagating through the pipe. The second meter portion includes a 
sensor for providing sound measurement signal indicative of the speed of sound propagating 
within the pipe. A processor provides a compensated meter measurement signal indicative 
of a measurement parameter corrected for entrained gas in the flow propagating through the 
pipe, in response to meter measurement signal and the sound measurement signal. 

The foregoing and other objects, features and advantages of the present invention 
will become more apparent in light of the following detailed description of exemplary 
embodiments thereof. 

Brief Description of the Drawings 

Fig. 1 is a schematic illustration of an apparatus having an array of sensors onto a 
pipe for measuring the volumetric flow and gas volume fraction of the mixture flowing in 
the pipe having entrained gas/air therein, in accordance with the present invention. 

Fig. 2 is a block diagram of an embodiment of the apparatus of Fig. 1, in accordance 
with the present invention. 

Fig. 3 is a functional flow diagram of an apparatus embodying the present invention 
that compensates the volumetric flow measurement of a volumetric flow meter, in 
accordance with the present invention. 
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Fig. 4 is a block diagram of an apparatus for measuring the speed of sound 
propagating through a process flow flowing within a pipe, in accordance with the present 
invention. 

Fig. 5 is a plot of Mixture Sound Speed as a function of gas volume fraction for a 
5 5% consistency slurry over a range of process pressures, in accordance with the present 
invention. 

Fig. 6 is a plot of Mixture Sound Speed a function of gas volume fraction for pure 
water and a 5% consistency slurry at 4 atm process pressure, in accordance with the present 
invention. 

1 0 Fig. 7 is a plot of Mixture Sound Speed as a function of gas volume fraction for 

different consistency slurry over a range of process pressures, in accordance with the 
present invention. 

Fig. 8 is a plot of Mixture Sound Speed a function of entrained air volume fraction 
for slurry at a process pressure, in accordance with the present invention. 
15 Fig. 9 is a K-w plot for acoustic field within 3 inch pipe containing -2% air by 

volume entrained in water flowing 240 gpm, in accordance with the present invention. 

Fig. 10 is a cross-sectional view of a pipe having a turbulent pipe flowing having 
coherent structures therein, in accordance with the present invention. 

Fig. 1 1 is a block diagram of an apparatus for measuring the vortical field of a 
20 process flow within a pipe, in accordance with the present invention. 

Fig. 12 a ko) plot of data processed from an apparatus embodying the present 
invention that illustrates slope of the convective ridge, and a plot of the optimization 
function of the convective ridge, in accordance with the present invention. 

Fig. 13 is a functional flow diagram of an apparatus embodying the present 
25 invention that compensates the volumetric flow measurement of an electromagnetic flow 
meter, in accordance with the present invention. 

Fig. 14 is a functional flow diagram of an apparatus embodying the present 
invention that compensates the consistency measurement of a consistency meter, in 
accordance with the present invention. 
30 Figs. 15 - 17 are configurations for an apparatus in accordance with the present 

invention. 
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Figs. 18 - 20 are plots of the output of an apparatus embodying the present invention 
for compensating a microwave consistency meter, in accordance with the present invention. 

Best Mode for Carrying Out the Invention 

5 Referring to Fig. 1, an apparatus, generally shown as 10, is provided to measure 

volumetric flow rate and gas volume fraction in liquids and mixtures (e.g. paper and pulp 
slurries or other solid liquid mixtures) having entrained gas therein (including air). The 
apparatus 10 in accordance with the present invention determines the speed at which sound 
propagates within a pipe 14 to measure entrained gas in liquids and/or mixtures 12. To 

10 simplify the explanation of the present invention the flow propagating through the pipe will 
be referred to as a mixture or slurry with the understanding that the flow may be a liquid or 
any other mixture having entrained gas therein. 

The following approach may be used with any technique that measures the sound 
speed of a fluid. However, it is particularly synergistic with sonar based volumetric flow 

15 meters such as described in U.S Patent Application, Serial No. (Cidra's Docket No. CC- 
0122A) and U.S. Patent Application, Serial No. 09/729,994 (Cidra's Docket No. CC-0297), 
filed December 4, 200, now US6,609,069, which are incorporated herein by reference, in 
that the sound speed measurement, and thus gas volume fraction measurement, can be 
accomplished using the same hardware as that required for the volumetric flow 

20 measurement. It should be noted, however, that the gas volume fraction measurement could 
be performed independently of a volumetric flow measurement, and would have utility as 
an important process measurement in isolation or in conjunction with other process 
measurements, which will be described in greater detail hereinafter. 

Fig. 2 is a block diagram 1 of the apparatus 10 of Fig. 1 that includes a device 2 for 

25 measuring the speed of sound (SOS) propagating within a pipe 14 and a device 3 for 
measuring the velocity of the mixture 12 within the pipe 14. A pressure sensor and/or 
temperature sensor 4 measures the pressure and/or temperature of the mixture flowing 
through the pipe. Alternatively, the pressure and/or temperature may be estimated rather 
than actually measured. In response to the speed of sound signal 5, the velocity 7 of the 

30 flow 12 and characteristics 6 of the flow (e.g., pressure and temperature), a processor 24 
determines the gas volume fraction (GVF) of the flow 12, the uncompensated volumetric 
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flow 9 of the mixture, and the volumetric flow 1 1 of the flow compensated for the entrained 
air therein. 

A flow chart 13 shown in Fig. 3 illustrates the function of the processor 24. As 
shown in Fig. 2, the inputs to the processor includes the speed of sound (SOS) within the 
5 pipe 14, the velocity of the mixture 12, and the pressure and temperature of the mixture. 
The fluid properties of the mixture (e.g., SOS and density) are determined knowing the 
pressure and temperature of the mixture. The gas volume fraction of the mixture (GVF) is 
determined using the SOS measurement and fluid properties, which will be described in 
greater detail hereinafter. The volumetric flow rate of the mixture (including the entrained 
10 gas) is determined using the velocity and knowing the cross-sectional area of the inner 
diameter of the pipe. The processor 24 provides a compensated volumetric flow 
measurement of the mixture by correcting the uncompensated volumetric flow rate using 
the void fraction of the air. For example, correction for void fraction of gas may be as 
follows, for a no slip, homogeneous flow model: 
1 5 Qair + Qliquid = Qmix 

Qair = GVFair * Q mix 

Qliquid = (1 -GVFair) Qmix 

Other models and corrections may be used to correct for gas volume fraction. 
Other information relating to the gas volume fraction in a fluid and the speed of 
20 sound (or sonic velocity) in the fluid, is described in "Fluid Mechanics and Measurements 
in two-phase flow Systems", Institution of mechanical engineers, proceedings 1969-1970 
Vol. 184 part 3C, Sept. 24-25 1969, Birdcage Walk, Westminster, London S.W. 1, England, 
which is incorporated herein by reference. 

Fig. 1 illustrates a schematic drawing of an embodiment of the present invention. 
25 The apparatus 10 includes a sensing device 16 comprising an array of pressure sensors (or 
transducers) 18-21 spaced axially along the outer surface 22 of a pipe 14, having a process 
flow propagating therein. The pressure sensors measure the unsteady pressures produced 
by acoustical and vortical disturbances within the pipe, which are indicative of the SOS 
propagating through the pipe and the velocity of the mixture 12. The output signals (P| - 
30 P N ) of the pressure sensors 1 8-2 1 are provided to the processor 24, which processes the 



-7- 

CiDRA Corporation Proprietary 

CC-0679PR EV US 



pressure measurement data and determines gas volume fraction (GVF), the uncompensated 
volumetric flow rate and the compensated volumetric flow rate, as described hereinbefore. 

In an embodiment of the present invention shown in Fig. 1, the apparatus 10 has at 
least four pressure sensors 18-21 disposed axially along the pipe 14 for measuring the 
5 unsteady pressure P r P N of the mixture 12 flowing therethrough. . Both measurement are 
derive by interpreting the unsteady pressure field within the process piping using multiple 
transducers displaced axially over ~ 2 diameters in length. The flow measurements can be 
performed using ported pressure transducers or clamp-on, strain-based sensors. 

The apparatus 10 has the ability to measure the gas volume fraction and volumetric 
10 flow rate using one or both of the following techniques described herein below: 

1) Determining the speed of sound of acoustical disturbances or sound waves 
propagating through the flow 12 using the array of pressure sensors 18-21, 
and/or 

2) Determining the velocity of vortical disturbances or "eddies" propagating 
1 5 through the flow 1 2 using the array of pressure sensors 18-21. 

Generally, the first technique measures unsteady pressures created by acoustical 
disturbances propagating through the flow 12 to determine the speed of sound (SOS) 
propagating through the flow. Knowing the pressure and/or temperature of the flow and the 
speed of sound of the acoustical disturbances, the processing unit 24 can determine the gas 

20 volume fraction of the mixture, as described and shown in Fig 3. 

The apparatus in Fig. 1 also contemplates providing one or more acoustic sources 27 
to enable the measurement of the speed of sound propagating through the flow for instances 
of acoustically quiet flow. The acoustic sources may be disposed at the input end of output 
end of the array of sensors 18-21, or at both ends as shown. One should appreciate that in 

25 most instances the acoustics sources are not necessary and the apparatus passively detects 
the acoustic ridge provided in the flow 12. The passive noise includes noise generated by 
pumps, valves, motors, and the turbulent mixture itself. 

The second technique measures the velocities associated with unsteady flow fields 
and/or pressure disturbances created by vortical disturbances or "eddies" 1 18 to determine 

30 the velocity of the flow 12. The pressure sensors 18-21 measure the unsteady pressures P r 
P N created by the vortical disturbances as these disturbances convect within the flow 12 
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through the pipe 14 in a known manner, as shown in Fig. 10. Therefore, the velocity of 
these vortical disturbances is related to the velocity of the mixture and hence the volumetric 
flow rate may be determined, as will be described in greater detail hereinafter. 

In one embodiment of the present invention as shown in Fig. 1, each of the pressure 
5 sensors 1 8-2 1 may include a piezoelectric film sensor to measure the unsteady pressures of 
the mixture 12 using either technique described hereinbefore. 

The piezoelectric film sensors include a piezoelectric material or film to generate an 
electrical signal proportional to the degree that the material is mechanically deformed or 
stressed. The piezoelectric sensing element is typically conformed to allow complete or 
10 nearly complete circumferential measurement of induced strain to provide a 

circumferential-averaged pressure signal. The sensors can be formed from PVDF films, co- 
polymer films, or flexible PZT sensors, similar to that described in "Piezo Film Sensors 
Technical Manual" provided by Measurement Specialties, Inc., which is incorporated herein 
by reference. A piezoelectric film sensor that may be used for the present invention is part 
1 5 number 1 - 1 002405-0, LDT4-028K, manufactured by Measurement Specialties, Inc. 

Piezoelectric film ("piezofilm"), like piezoelectric material, is a dynamic material 
that develops an electrical charge proportional to a change in mechanical stress. 
Consequently, the piezoelectric material measures the strain induced within the pipe 14 due 
to unsteady pressure variations (e.g., vortical and/or acoustical) within the process mixture 
20 12. Strain within the pipe is transduced to an output voltage or current by the attached 

piezoelectric sensor. The piezoelectrical material or film may be formed of a polymer, such 
as polarized fluoropolymer, polyvinylidene fluoride (PVDF). The piezoelectric film 
sensors are similar to that described in U.S. Patent Application Serial No. (CiDRA Docket 
No. CC-0676), which is incorporated herein by reference. 
25 The apparatus 10 of the present invention may be configured and programmed to 

measure and process the detected unsteady pressures Pi(t) - P N (t) created by acoustic waves 
and/or vortical disturbances, respectively, propagating through the mixture to determine the 
SOS within the pipe 14 and the velocity of the mixture 12. One such apparatus 1 10 is 
shown in Fig. 4 that measures the speed of sound (SOS) of one-dimensional sound waves 
30 propagating through the mixture to determine the gas volume fraction o f the mixture. It is 
known that sound propagates through various mediums at various speeds in such fields as 
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SONAR and RADAR fields. The speed of sound propagating through the pipe and mixture 
12 may be determined using a number of known techniques, such as those set forth in U.S. 
Patent Application Serial No. 09/344,094, entitled "Fluid Parameter Measurement in Pipes 
Using Acoustic Pressures", filed June 25, 1999, now US 6,354,147; U.S. Patent Application 
5 Serial No. 09/729,994, filed December 4, 2002, now US 6,609,069; U.S. Patent Application 
Serial No. 09/997,221, filed November 28, 2001, now US 6,587,798; and U.S. Patent 
Application Serial No. 10/007,749, entitled "Fluid Parameter Measurement in Pipes Using 
Acoustic Pressures", filed November 7, 2001, each of which are incorporated herein by 
reference. 

10 In accordance with the present invention, the speed of sound propagating through 

the mixture 12 is measured by passively listening to the flow with an array of unsteady 
pressure sensors to determine the speed at which one-dimensional compression waves 
propagate through the mixture 12 contained within the pipe 14. 

As shown in Fig. 4, the apparatus 1 10 has an array of at least three acoustic pressure 

15 sensors 1 15,1 16,1 17, located at three locations Xi,x 2 ,x 3 axially along the pipe 14. One will 
appreciate that the sensor array may include more than three pressure sensors as depicted by 
pressure sensor 1 18 at location x N . The pressure generated by the acoustic waves may be 
measured through pressure sensors 115-118. The pressure sensors 15-18 provide 
pressure time-varying signals Pi(t),P 2 (t),P 3 (t),P N (t) on lines 120,121,122,123 to a signal 

20 processing unit 130 to known Fast Fourier Transform (FFT) logics 126,127,128,129, 

respectively. The FFT logics 126- 129 calculate the Fourier transform of the time-based 
input signals Pi(t) - P N (t) and provide complex frequency domain (or frequency based) 
signals Pi(cd),P 2 (co),P3(o>),Pn(o)) on lines 132,133,134,135 indicative of the frequency 
content of the input signals. Instead of FFT's, any other technique for obtaining the 

25 frequency domain characteristics of the signals Pi(t) - P N (t), may be used. For example, the 
cross-spectral density and the power spectral density may be used to form a frequency 
domain transfer functions (or frequency response or ratios) discussed hereinafter. 

The frequency signals P^co) - P N (co) are fed to a mix -Mx Calculation Logic 138 which 
provides a signal to line 40 indicative of the speed of sound of the mixture a™* (discussed 

30 more hereinafter). The a m j X signal is provided to map (or equation) logic 142, which 
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converts to a percent composition of a mixture and provides a %Comp signal to line 44 
indicative thereof (as discussed hereinafter). 

More specifically, for planar one-dimensional acoustic waves in a homogenous 
mixture, it is known that the acoustic pressure field P(x,t) at a location x along the pipe 14, 
5 where the wavelength X of the acoustic waves to be measured is long compared to the 
diameter d of the pipe 14 (i.e., X/d »1), may be expressed as a superposition of a right 
traveling wave and a left traveling wave, as follows: 

P(x, t) = (Ae -f Be + ik ' x }> im Eq. 1 

where A,B are the frequency-based complex amplitudes of the right and left traveling 
10 waves, respectively, x is the pressure measurement location along a pipe 14, <o is frequency 
(in rad/sec, where G)=2rcf), and kr,ki are wave numbers for the right and left traveling waves, 
respectively, which are defined as: 



\ + M x 



and hi s 



O) 



\ a mixJ 



\-M x 



Eq.2 



where a m j X is the speed of sound of the mixture in the pipe, co is frequency (in rad/sec), and 
15 M x is the axial Mach number of the flow of the mixture within the pipe, where: 

M x mZs*L Eq<3 

a mix 

where Vmix is the axial velocity of the mixture. For non-homogenous mixtures, the axial 
Mach number represents the average velocity of the mixture and the low frequency acoustic 
field description remains substantially unaltered. 

20 The data from the array of sensors 115-118 may be processed in any domain, 

including the frequency/spatial domain, the temporal/spatial domain, the temporal/wave- 
number domain or the wave-number/frequency (k-©) domain. As such, any known array 
processing technique in any of these or other related domains may be used if desired, 
similar to the techniques used in the fields of SONAR and RADAR. 

25 Also, some or all of the functions within the processor 130 may be implemented in 

software (using a microprocessor or computer) and/or firmware, or may be implemented 
using analog and/or digital hardware, having sufficient memory, interfaces, and capacity to 
perform the functions described herein. 
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Acoustic pressure sensors 115-118 sense acoustic pressure signals that, as 
measured, are lower frequency (and longer wavelength) signals than those used for 
ultrasonic probes of the prior art, and thus the current invention is more tolerant to 
inhomogeneities in the flow, such as time and space domain inhomogeneities within the 



be known or arbitrary and that as few as two sensors are required if certain information is 
known about the acoustic properties of the process flow 12. The pressure sensors are 
spaced sufficiently such that the entire length of the array (aperture) is at least a significant 
fraction of the measured wavelength of the acoustic waves being measured. The acoustic 
wavelength is a function of the type or characteristics of flow 12. 

Based on the above discussion, one may use a short length scale aperture to measure 
the sound speed. 

The characteristic acoustic length scale is: X=c/f; where c is the speed of sound in a 
mixture, f is frequency and X is wavelength. 

If Aperture = L and if L/X is approx. constant. 

Then Lwater/Xwater = Lwater * f7C wa ter « Lgvf *ffc G vF 

Therefore: L G vf = Lwater (C G vF/C W ater); where GVF is gas volume fraction. 

Thus for SOS of water (Cwater = 5,000 ft/sec) , and SOS of the Gas volume 
fraction (C GVF = 500 ft/sec) and a length aperture of L water = 5ft (which we have shown 
is sufficient to accurately measure the SOS of water), the length aperture for a gas volume 
fraction Lgvf would be about 0.5 feet. 

The sound speed of a mixture can be related to volumetric phase fraction (<|>j) of the 
components and the sound speed (a) and densities (p) of the component through the Wood 
equation. 



One dimensional compression waves propagating within a mixture 12 contained 
within a pipe 14 exert an unsteady internal pressure loading on the pipe. The degree to 



flow. 



It is within the scope of the present invention that the pressure sensor spacing may 
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which the pipe displaces as a result of the unsteady pressure loading influences the speed of 
propagation of the compression wave. The relationship among the infinite domain speed of 
sound and density of a mixture; the elastic modulus (E), thickness (t), and radius (R) of a 
vacuum-backed cylindrical conduit; and the effective propagation velocity (a^) for one 
5 dimensional compression is given by the following expression: 




Note: "vacuum backed" as used herein refers to a situation in which the fluid 

10 surrounding the pipe externally has negligible acoustic impedance compared to that of the 
mixture internal to the pipe 14. For example, meter containing a typical water and pulp 
slurry immersed in air at standard atmospheric conditions satisfies this condition and can be 
considered "vacuum-backed". 

The mixing rule essentially states that the compressibility of a mixture (l/(p a 2 )) is 

15 the volumetrically-weighted average of the compressibilities of the components. For gas / 
liquid mixtures 12 at pressure and temperatures typical of paper and pulp industry, the 
compressibility of gas phase is orders of magnitudes greater than that of the liquid. Thus, 
the compressibility of the gas phase and the density of the liquid phase primarily determine 
mixture sound speed., and as such, it is necessary to have a good estimate of process 

20 pressure to interpret mixture sound speed in terms of volumetric fraction of entrained air. 
The effect of process pressure on the relationship between sound speed and entrained air 
volume fraction is shown in Fig. 5. 

Conversely, however, detailed knowledge of the liquid / slurry is not required for 
entrained air measurement. Variations in liquid density and compressibility with changes in 

25 consistency have a negligible effect on mixture sound speed compared to the presence of 
entrained air. Fig. 6 shows the mixture sound speed as a function of entrained air volume 
fraction for two slurries, one with 0% wood fiber and the other with 5% wood fiber by 
volume. As shown, the relationship between mixture sound speed and gas volume fraction 
is essentially indistinguishable for the two slurries. Furthermore, mixture sound speed is 
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shown to an excellent indicator of gas volume fraction, especially for the trace to moderate 
amounts of entrained air, from 0 to 5% by volume, typically encountered in the paper and 
pulp industry. 

For paper and pulp slurries, the conditions are such that for slurries with non- 
5 negligible amounts of entrained gas, say < 0.01%, the compliance of standard industrial 
piping (Schedule 10 or 40 steel pipe) is typically negligible compared to that of the 
entrained air. 

Figs. 7 and 8 above show the relationship between sound speed and entrained air for 
slurries 12 with pulp contents representative of the range used in the paper and pulp 

10 industry. Referring to Fig. 7, two slurry consistencies are shown; representing the lower 
limit, a pure water mixture is considered, and representing the higher end of consistencies, a 
5% pulp / 95% water slurry is considered. Since the effect of entrained air on the sound 
speed of the mixture is highly sensitive to the compressibility of the entrained air, the effect 
of the entrained air is examined at two pressures, one at ambient representing the lower 

15 limit of pressure, and one at 4 atmospheres representing a typical line pressure in a paper 
process. As shown, the consistency of the liquid slurry 12, i.e., the pulp content, has little 
effect on the relationship between entrained air volume fraction and mixture sound speed. 
This indicates that an entrained air measurement could be accurately performed, within 
0.01% or so, with little or no knowledge of the consistency of the slurry. The chart does 

20 show a strong dependence on line pressure. Physically, this effect is linked to the 

compressibility of the air, and thus, this indicates that reasonable estimates of line pressure 
and temperature would be required to accurately interpret mixture sound speed in terms of 
entrained air gas volume fraction. 

Fig. 7 also shows that for the region of interest, from roughly 1% entrained air to 

25 roughly 5% entrained air, mixture sound speeds (amix) are quite low compare to the liquid- 
only sound speeds. In this example, the sound speed of the pure water and the 5% pulp 
slurry were calculated, based on reasonable estimates of the constituent densities and 
compressibilities, to be 1524 m/s and 1541 m/s, respectively. The sound speed of these 
mixtures with 1% to 5% entrained air at typical operating pressure (latm to 4 atms) are on 

30 the order of 100 m/sec. The implication of these low sound speed is that the mixture sound 
speed could be accurately determined with a array of sensors, ie using the methodology 
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described in aforementioned U.S Patent Applications, Serial No. (Cidra's Docket No. CC- 
0066A), and/or Serial No. (Cidra's Docket No. CC-0066B), with an aperture that is similar, 
or identical, to an array of sensors that would be suitable to determine the convection 
velocity, using the methodology described in aforementioned U.S Patent Application, Serial 
5 No. (Cidra's Docket No. CC-0122A), which is incorporated herein by reference. Thus, 
performing a volumetric flow measurement and an entrained air volumetric flow 
measurement using the convection velocity and mixture sound speed simultaneously, with 
the same sensor array would provide functionality currently unavailable to the paper and 
pulp industry. 

10 As mentioned earlier, the relationship between mixture sound speed and entrained 

air in bubbly liquids 12 is well established. However, as will be developed below, in 
bubbly flows, these relations are only applicable for the propagation of relatively low 
frequency, long wavelength sound. While this restriction does not present any significant 
obstacles for the sonar meter, it does present significant challenges to ultrasonic sound 

1 5 speed measurement devices. 

Ultrasonic meters typically operate in 100 KHz to several MHz frequency range. 
For these meters, entrained air bubbles have length scales on the same order as the acoustic 
waves generated by the ultrasonic meters. The posed several problems. Firstly, the bubbles 
scatter the ultrasonic waves, impairing the ability of the ultrasonic meter to perform a sound 

20 speed measurement. Also, ultrasonic meters rely on information derived from only a small 
fraction of the cross sectional area of the pipe to be representative of the entire cross 
section, an assumption that breaks down for flows with inhogenieties on the same length 
scale as the ultrasonic wavelength. 

For the sound speed measurement, the apparatus 1 10 utilizes similar processing 

25 algorithms as those employed for the volumetric flow measurement. As with convective 
disturbances (which is described in greater detail hereinafter), the temporal and spatial 
frequency content of sound propagating within the process piping is related through a 
dispersion relationship. 
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As before, k is the wave number, defined as k=27i/A., © is the temporal frequency in 
rad/sec, and a,™ is the speed at which sound propagates within the process piping. Unlike 
disturbances, which convect with the flow, however, sound generally propagates in both 
directions, with and against the mean flow. For these cases, the acoustic power is located 
5 along two acoustic ridges, one for the sound traveling with the flow at a speed of a mix + V mix 
and one for the sound traveling against the flow at a speed of a mix - V mix . 

Fig. 9 shows a k-o> plot generated for acoustic sound field recorded from water 
flowing at a rate of 240 gpm containing -2% entrained air by volume in a 3 in, schedule 10, 
stainless steel pipe. The k-<o plot was constructed using data from an array of strain-based 

10 sensors attached to the outside of the pipe. Two acoustic ridges are clearly evident. Based 
on the slopes of the acoustic ridges, the sound speed for this for this mixture was 330 ft/sec 
( 1 OOm/s), consistent with that predicted by the Wood equation. Note that adding 2% air by 
volume reduces the sound speed of the bubbly mixture to less than 10% of the the sound 
speed of single phase water. 

15 While the sonar-based flow meter using an array of sensors to measure the speed of 

sound of an acoustic wave propagating through the mixture, one will appreciate that any 
means for measuring the speed of sound of the acoustic wave may used to determine the 
entrained air volume fraction of the mixture/fluid. 

The apparatus 1 10 further includes the ability to measure of volumetric flow rate of 

20 the mixture by comparing the difference of the speed of one dimensional sound waves 
propagating with and against the mean flow. 

This method of determining the volumetric flow rate of the flow 12 relies on the 
interaction of the mean flow with the acoustic pressure field. The interaction results in 
sound waves propagating with the mean flow traveling at the speed of sound (if the 

25 vapor/liquid mixture were not flowing) plus the convection velocity and, conversely, sound 
waves traveling against the mean flow propagating at the speed of sound minus the 
convection velocity. That is, 

= a mix + u 
30 a L = a mix - u 
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where a R = velocity of a right traveling acoustic wave relative to a stationary observer (i.e. 
the tube 14), a L = velocity of a left traveling acoustic wave apparent to a stationary 
observer, a m i x = speed of sound traveling through the mixture (if the mixture was not 
flowing) and u = the mean flow velocity (assumed to be flowing from left to right in this 
5 instance). Combining these two equations yields an equation for the mean velocity, 




Therefore, by measuring the propagation velocity of acoustic waves in both directions 

10 relative to the pipe 14 as described hereinbefore, the mean flow velocity can be calculated 
by multiplying the mean flow velocity by the cross-sectional area of the pipe 14. 

Further, Fig. 9 illustrates the ability of the present invention to determine the 
velocity of a fluid moving in a pipe. The color contours represent the relative signal power 
at all combinations of frequency and wavenumber. The highest power "ridges" represent 

15 the acoustic wave with slope of the ridges equal to the propagation speed. The dashed lines 
show the best-fit two-variable maximization of the power with the two variables being 
sound speed and flow velocity. The right-side ridge represents the acoustic wave traveling 
in the same direction as the bulk flow and therefore its slope is steeper than the left-side 
ridge that represents the acoustic wave traveling in the opposite direction of the flow. This 

20 indicates that the acoustic wave traveling in the same direction of the flow is traveling faster 
than the acoustic wave traveling in the opposite direction of the flow relative to the 
stationary sensors located on the probe. 

As discussed hereinbefore, the apparatus 10 of Fig. 1 embodying the present 
invention also includes the ability to measure volumetric flow rate of the mixture by 

25 measuring the unsteady pressures generated by vortical disturbances 1 88 propagating in the 
mixture 12 (see Fig. 10). The apparatus 10 uses one or both of the following techniques to 
determine the convection velocity of the vortical disturbances within the process flow 12 
by: 

1) Cross-correlating unsteady pressure variations using an array of unsteady pressure 
30 sensors. 
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2) Characterizing the convective ridge of the vortical disturbances using an array of 
unsteady pressure sensors. 



5 To measure volumetric flow, the sonar meter characterizes speed at which coherent 

vortical structures convect past an axial array of sensors using beam forming techniques 
developed over several decades for underwater acoustic application. Coherent structures 
are an inherent feature of turbulent boundary layers present in all turbulent flows. Unlike 
conventional vortex shedding meters, no internal geometry is required to generate these 
10 structures. 

The overwhelming majority of industrial process flows involve turbulent flow 12. 
Turbulent fluctuations within the process flow govern many of the flow properties of 
practical interest including the pressure drop, heat transfer, and mixing. For engineering 
applications, considering only the time-averaged properties of turbulent flows is often 

1 5 sufficient for design purposes. For sonar based array processing flow metering technology, 
understanding the time-averaged velocity profile in turbulent flow 12 provides a means to 
interpret the relationship between speed at which coherent structures 118 convect and the 
volumetrically averaged flow rate. 

Turbulent pipe flows 12 are highly complex flows. Predicting the details of any 

20 turbulent flow is problematic, however, much is known regarding the statistical properties 
of the flow. For instance, turbulent flows contain self-generating, coherent vortical 
structures often termed "turbulent eddies". The maximum length scale of these eddies is set 
by the diameter of the pipe 14. These structures remain coherent for several tube diameters 
downstream, eventually breaking down into progressively smaller eddies until the energy is 

25 dissipated by viscous effects. 

Experimental investigations have established that eddies generated within turbulent 
boundary layers convect at roughly 80% of maximum flow velocity. For tube flows, this 
implies that turbulent eddies will convect at approximately the volumetrically averaged 
flow velocity within the pipe 14. The precise relationship between the convection speed of 

30 turbulent eddies and the flow rate for each class of meters can be calibrated empirically as 
described below. 
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Fig. 10 illustrates the relevant flow features of turbulent pipe flow 12 along with a 
axial array of sensors 18-21. As shown, the time-averaged axial velocity is a function of 
radial position, from zero the wall to a maximum at the centerline of the pipe. The flow 12 
near the wall is characterized by steep velocity gradients and transitions to relatively 
uniform core flow near the center of the pipe 14. Vortical structures, often termed turbulent 
eddies, are superimposed over time averaged velocity profile. These coherent structures 
contain temporally and spatially random fluctuations with magnitudes typically less than 
10% percent of the mean flow velocity and are carried along with the mean flow. 
Experimental investigations have established that eddies generated within turbulent 
boundary layers remain coherent for several pipe diameters and convect at roughly 80% of 
maximum flow velocity (Schlichting, 1979). 

From a volumetric flow measurement perspective, the volumetrically averaged flow 
velocity is of interest. The volumetrically averaged flow velocity, defined as the total 
volumetric flow rate, Q, divided by the cross sectional area of the conduit, A, is a useful, but 
arbitrarily defined property of the flow. In fact, given the velocity profile within the pipe, 
little flow is actually moving at this speed. The precise relationship between the convection 
speed of turbulent eddies and the flow rate is determined experimentally through calibration 
for each. 

The Reynolds number (Re), based on pipe diameter (D), characterizes many of the 
engineering properties of the flow. The Reynolds number is a non-dimensional ratio 
representing the relative importance of inertial forces to viscous forces within a flow: 



inertial ^ fa VD 

Re = forces = ~- = 

viscous d u v 

Where p is the fluid density, p is the dynamic viscosity, U is the volumetrically averaged 
flow velocity and v (= p /p) is the kinematic viscosity. 

The critical Reynolds number for pipe flows, above which flows are considered 
turbulent, is -2300. Most flows in the paper and pulp industry have Reynolds number 
ranging from one hundred thousand to several million, well within the turbulent regime. In 
addition to demarcating a boundary between laminar and turbulent flow regimes, the 
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Reynolds number is a similarity parameter for pipe flows, i.e. flows in geometrically similar 
pipes with the same Reynolds number are dynamically similar (Schlichting p. 12). 

The apparatus 170 of Fig. 1 1 determines the convection velocity of the vortical 
disturbances within the flow by cross correlating unsteady pressure variations using an array 
of unsteady pressure sensors, similar to that shown in U.S. Patent Application Serial No. 
10/007,736, filed November 8, 2001, entitled "Flow Rate Measurement Using Unsteady 
Pressures", which is incorporated herein by reference. 

Referring to Fig. 1 1, the apparatus 170 includes a sensing section 172 along a pipe 
14 and a signal processing unit 174. The pipe 14 has two measurement regions 176,178 
located a distance AX apart along the pipe 14. At the first measurement region 176 are two 
unsteady (or dynamic or ac) pressure sensors 180,182, located a distance Xi apart, capable 
of measuring the unsteady pressure in the pipe 14, and at the second measurement region 
178, are two other unsteady pressure sensors 84,86, located a distance X 2 apart, capable of 
measuring the unsteady pressure in the pipe 14. Each pair of pressure sensors 180,182 and 
184,186 act as spatial filters to remove certain acoustic signals from the unsteady pressure 
signals, and the distances X!,X 2 are determined by the desired filtering characteristic for 
each spatial filter, as discussed hereinafter. 

The apparatus 170 of the present invention measures velocities associated with 
unsteady flow fields and/or pressure disturbances represented by 188 associated therewith 
relating to turbulent eddies (or vortical flow fields), inhomogeneities in the flow, or any 
other properties of the flow, liquid, vapor, or pressure, having time varying or stochastic 
properties that are manifested at least in part in the form of unsteady pressures. The vortical 
flow fields are generated within the flow of the pipe 14 by a variety of non-discrete sources 
such as remote machinery, pumps, valves, elbows, as well as the fluid or mixture flow itself. 
It is this last source, the fluid flowing within the pipe, that is a generic source of vortical 
flow fields primarily caused by the shear forces between the flow 12 and the wall of the 
tube that assures a minimum level of disturbances for which the present invention takes 
unique advantage. The flow generated vortical flow fields generally increase with mean 
flow velocity and do not occur at any predeterminable frequency. As such, no external 
discrete vortex-generating source is required within the present invention and thus may 
operate using passive detection. It is within the scope of the present that the pressure sensor 
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spacing may be known or arbitrary and that as few as two sensors are required if certain 
information is known about the acoustic properties of the system as will be more fully 
described herein below. 

The vortical flow fields 188 are, in general, comprised of pressure disturbances 
5 having a wide variation in length scales and which have a variety of coherence length scales 
such as that described in the reference "Sound and Sources of Sound", A. P.Dowling et al, 
Halsted Press, 1983, which is incorporated by reference to the extend of understanding the 
invention. Certain of these vortical flow fields 188 convect at or near, or related to the 
mean velocity of at least one of the elements within a mixture flowing through the pipe 14. 

10 The vortical pressure disturbances 1 88 that contain information regarding convection 
velocity have temporal and spatial length scales as well as coherence length scales that 
differ from other disturbances in the flow. The present invention utilizes these properties to 
preferentially select disturbances of a desired axial length scale and coherence length scale 
as will be more fully described hereinafter. For illustrative purposes, the terms vortical flow 

15 field and vortical pressure field will be used to describe the above-described group of 

unsteady pressure fields having temporal and spatial length and coherence scales described 
herein. 

Also, some or all of the functions within the signal processing unit 174 may be 
implemented in software (using a microprocessor or computer) and/or firmware, or may be 
20 implemented using analog and/or digital hardware, having sufficient memory, interfaces, 
and capacity to perform the functions described herein. 

In particular, in the processing unit 174, the pressure signal Pi(t) on the line 190 is 
provided to a positive input of a summer 200 and the pressure signal P 2 (t) on the line 191 is 
provided to a negative input of the summer 200. The output of the summer 200 is provided 
25 to line 204 indicative of the difference between the two pressure signals Pi,P 2 (e.g., P r 

P 2 =Pasl). 

The pressure sensors 180,182 together with the summer 200 create a spatial filter 
176. The line 204 is fed to bandpass filter 208, which passes a predetermined passband of 
frequencies and attenuates frequencies outside the passband. In accordance with the present 
30 invention, the passband of the filter 208 is set to filter out (or attenuate) the dc portion and 
the high frequency portion of the input signals and to pass the frequencies therebetween. 
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Other passbands may be used in other embodiments, if desired. Passband filter 208 
provides a filtered signal Pasfl on a line 212 to Cross-Correlation Logic 216, described 
hereinafter. 

The pressure signal P 3 (t) on the line 192 is provided to a positive input of a summer 
5 202 and the pressure signal P 4 (t) on the line 193 is provided to a negative input of the 
summer 202. The pressure sensors 83,84 together with the summer 202 create a spatial 
filter 178. The output of the summer 202 is provided on a line 206 indicative of the 
difference between the two pressure signals P 3 ,P 4 (e.g., P 3 -P 4 = Pas2). The line 206 is fed to a 
bandpass filter 210, similar to the bandpass filter 108 discussed hereinbefore, which passes 

10 frequencies within the passband and attenuates frequencies outside the passband. The filter 
210 provides a filtered signal Pasf2 on a line 214 to the Cross-Correlation Logic 216. The 
signs on the summers 200,202 may be swapped if desired, provided the signs of both 
summers are swapped together. In addition, the pressure signals Pi,P 2 ,P3,P4 may be scaled 
prior to presentation to the summers 200,202. 

15 The Cross-Correlation Logic 216 calculates a known time domain cross-correlation 

between the signals Pasn and P^o on the lines 212,214, respectively, and provides an output 
signal on a line 218 indicative of the time delay x it takes for an vortical flow field 188 (or 
vortex, stochastic, or vortical structure, field, disturbance or perturbation within the flow) to 
propagate from one sensing region 176 to the other sensing region 178. Such vortical flow 

20 disturbances, as is known, are coherent dynamic conditions that can occur in the flow which 
substantially decay (by a predetermined amount) over a predetermined distance (or 
coherence length) and convect (or flow) at or near the average velocity of the fluid flow. As 
described above, the vortical flow field 188 also has a stochastic or vortical pressure 
disturbance associated with it. In general, the vortical flow disturbances 188 are distributed 

25 throughout the flow, particularly in high shear regions, such as boundary layers (e.g., along 
the inner wall of the tube 14) and are shown herein as discrete vortical flow fields 188. 
Because the vortical flow fields (and the associated pressure disturbance) convect at or near 
the mean flow velocity, the propagation time delay x is related to the velocity of the flow by 
the distance AX between the measurement regions 176,178, as discussed hereinafter. 

30 Referring to Fig. 15, a spacing signal AX on a line 220 indicative of the distance AX 

between the sensing regions 176,178 is divided by the time delay signal x on the line 218 by 
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a divider 222 which provides an output signal on the line 196 indicative of the convection 
velocity U c (t) of the saturated vapor/liquid mixture flowing in the pipe 14, which is related 
to (or proportional to or approximately equal to) the average (or mean) flow velocity Uf(t) 
of the flow 12, as defined below: 

5 U c (t) = AX/xocUf(t) Eq. 1 

The present invention uses temporal and spatial filtering to precondition the pressure 
signals to effectively filter out the acoustic pressure disturbances P aC oustic and other long 
wavelength (compared to the sensor spacing) pressure disturbances in the tube 14 at the two 

10 sensing regions 1 76, 1 78 and retain a substantial portion of the vortical pressure disturbances 
Pvorticai associated with the vortical flow field 188 and any other short wavelength (compared 
to the sensor spacing) low frequency pressure disturbances P ot her. In accordance with the 
present invention, if the low frequency pressure disturbances P ot her are small, they will not 
substantially impair the measurement accuracy of Pvorticai- 

1 5 The second technique of determining the convection velocity of the vortical 

disturbances within the flow 12 is by characterizing the convective ridge of the vortical 
disturbances using an array of unsteady pressure sensors, similar to that shown in U.S. 
Patent Application Serial No. 09/729,994, filed December 4, 2000, entitled "Method and 
Apparatus for Determining the Flow Velocity Within a Pipe", which is incorporated herein 

20 by reference. 

The sonar flow metering methodology uses the convection velocity of coherent 
structure with turbulent pipe flows 12 to determine the volumetric flow rate. The convection 
velocity of these eddies 188 is determined by applying sonar arraying processing techniques 
to determine the speed at which the eddies convect past an axial array of unsteady pressure 
25 measurements distributed along the pipe 14. 

The sonar-based algorithms determine the speed of the eddies 188 by characterizing 
both the temporal and spatially frequency characteristics of the flow field. For a series of 
coherent eddies convecting past a fixed array of sensors, the temporal and spatial frequency 
content of pressure fluctuations are related through the following relationship: 

30 
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03= — 
U. 



convect 



Here k is the wave number, defined as k=2n/X and has units of 1 /length, co is the temporal 
frequency in rad/sec, and Uconvect is the convection velocity. Thus, the shorter the 
wavelength (larger k) is, the higher the temporal frequency. 
5 In sonar array processing, the spatial / temporal frequency content of time stationary 

sound fields are often displayed using "k-co plots". K-co plots are essentially three- 
dimensional power spectra in which the power of a sound field is decomposed into bins 
corresponding to specific spatial wave numbers and temporal frequencies. On a k-co plot, 
the power associated with a pressure field convecting with the flow is distributed in regions, 

10 which satisfies the dispersion relationship developed above. This region is termed "the 
convective ridge" (Beranek, 1992) and the slope of this ridge on a k-w plot indicates the 
convective velocity of the pressure field. This suggests that the convective velocity of 
turbulent eddies, and hence flow rate within a tube, can be determined by constructing a k-co 
plot from the output of a phased array of sensor and identifying the slope of the convective 

15 ridge. 

Fig. 12 shows an example of a k-co plot generated from a phased array of pressure 
sensors. The power contours show a well-defined convective ridge. A parametric 
optimization method was used to determine the "best" line representing the slope of the 
convective ridge 200. For this case, a slope of 14.2 ft/sec was determined. The 

20 intermediate result of the optimization procedure is displayed in the insert, showing that 
optimized value is a unique and well-defined optima. 

The k-w plot shown in Fig. 12 illustrates the fundamental principle behind sonar 
based flow measure, namely that axial arrays of pressure sensors can be used in conjunction 
with sonar processing techniques to determine the speed at which naturally occurring 

25 turbulent eddies convect within a pipe. 

Another embodiment of the present invention include a pressure sensor such as pipe 
strain sensors, accelerometers, velocity sensors or displacement sensors, discussed 
hereinafter, that are mounted onto a strap to enable the pressure sensor to be clamped onto 
the pipe. The sensors may be removable or permanently attached via known mechanical 

30 techniques such as mechanical fastener, spring loaded, clamped, clam shell arrangement, 
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strapping or other equivalents. These certain types of pressure sensors, it may be desirable 
for the pipe 12 to exhibit a certain amount of pipe compliance. 

Instead of single point pressure sensors 18-21, at the axial locations along the pipe 
12, two or more pressure sensors may be used around the circumference of the pipe 12 at 
each of the axial locations. The signals from the pressure sensors around the circumference 
at a given axial location may be averaged to provide a cross-sectional (or circumference) 
averaged unsteady acoustic pressure measurement. Other numbers of acoustic pressure 
sensors and annular spacing may be used. Averaging multiple annular pressure sensors 
reduces noises from disturbances and pipe vibrations and other sources of noise not related 
to the one-dimensional acoustic pressure waves in the pipe 12, thereby creating a spatial 
array of pressure sensors to help characterize the one-dimensional sound field within the 
pipe 12. 

The pressure sensors 18-21 of Fig. 1 described herein may be any type of pressure 
sensor, capable of measuring the unsteady (or ac or dynamic ) pressures within a pipe 14, 
such as piezoelectric, optical, capacitive, resistive (e.g., Wheatstone bridge), accelerometers 
(or geophones), velocity measuring devices, displacement measuring devices, etc. If optical 
pressure sensors are used, the sensors 18-21 may be Bragg grating based pressure sensors, 
such as that described in US Patent Application, Serial No. 08/925,598, entitled " High 
Sensitivity Fiber Optic Pressure Sensor For Use In Harsh Environments", filed Sept. 8, 
1997, now U.S. Patent 6,016,702, and in US Patent Application, Serial No. 10/224,821, 
entitled " Non-Intrusive Fiber Optic Pressure Sensor for Measuring Unsteady Pressures 
within a Pipe", which are incorporated herein by reference. In an embodiment of the 
present invention that utilizes fiber optics as the pressure sensors 14 they may be connected 
individually or may be multiplexed along one or more optical fibers using wavelength 
division multiplexing (WDM), time division multiplexing (TDM), or any other optical 
multiplexing techniques. 

In certain embodiments of the present invention, a piezo-electronic pressure 
transducer may be used as one or more of the pressure sensors 15-18 and it may measure 
the unsteady (or dynamic or ac) pressure variations inside the tube 14 by measuring the 
pressure levels inside of the tube. These sensors may be ported within the pipe to make 
direct contact with the mixture 12. In an embodiment of the present invention, the sensors 
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14 comprise pressure sensors manufactured by PCB Piezotronics. In one pressure sensor 
there are integrated circuit piezoelectric voltage mode-type sensors that feature built-in 
microelectronic amplifiers, and convert the high-impedance charge into a low-impedance 
voltage output. Specifically, a Model 106B manufactured by PCB Piezotronics is used 
which is a high sensitivity, acceleration compensated integrated circuit piezoelectric quartz 
pressure sensor suitable for measuring low pressure acoustic phenomena in hydraulic and 
pneumatic systems. It has the unique capability to measure small pressure changes of less 
than 0.001 psi under high static conditions. The 106B has a 300 mV/psi sensitivity and a 
resolution of 91 dB (0.0001 psi). 

The pressure sensors incorporate a built-in MOSFET microelectronic amplifier to 
convert the high-impedance charge output into a low-impedance voltage signal. The sensor 
is powered from a constant-current source and can operate over long coaxial or ribbon cable 
without signal degradation. The low-impedance voltage signal is not affected by 
triboelectric cable noise or insulation resistance-degrading contaminants. Power to operate 
integrated circuit piezoelectric sensors generally takes the form of a low-cost, 24 to 27 
VDC, 2 to 20 mA constant-current supply. A data acquisition system of the present 
invention may incorporate constant-current power for directly powering integrated circuit 
piezoelectric sensors. 

Most piezoelectric pressure sensors are constructed with either compression mode 
quartz crystals preloaded in a rigid housing, or unconstrained tourmaline crystals. These 
designs give the sensors microsecond response times and resonant frequencies in the 
hundreds of kHz, with minimal overshoot or ringing. Small diaphragm diameters ensure 
spatial resolution of narrow shock waves. 

The output characteristic of piezoelectric pressure sensor systems is that of an AC- 
coupled system, where repetitive signals decay until there is an equal area above and below 
the original base line. As magnitude levels of the monitored event fluctuate, the output 
remains stabilized around the base line with the positive and negative areas of the curve 
remaining equal. 

It is also within the scope of the present invention that any strain sensing technique 
may be used to measure the variations in strain in the pipe, such as highly sensitive 
piezoelectric, electronic or electric, strain gages and piezo-resistive strain gages attached to 
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the pipe 12. Other strain gages include resistive foil type gages having a race track 
configuration similar to that disclosed U.S. Patent Application Serial No. 09/344,094, filed 
June 25, 1999, now US 6,354,147, which is incorporated herein by reference. The 
invention also contemplates strain gages being disposed about a predetermined portion of 
5 the circumference of pipe 12. The axial placement of and separation distance AX b AX 2 
between the strain sensors are determined as described herein above. 

It is also within the scope of the present invention that any other strain sensing 
technique may be used to measure the variations in strain in the tube, such as highly 
sensitive piezoelectric, electronic or electric, strain gages attached to or embedded in the 
10 tube 14. 

Note that this entrained air or gas volume fraction measurement GVFair may be 
used with any flow meter to correct for errors introduced into a measurement by entrained 
air. For instance, an electromagnetic flow meter will show an error when entrained air exists 

15 in the mixture. The present invention may be used to correct for this error. 

Referring to 13, the present invention contemplates an apparatus 300 for providing a 
volumetric flow measurement 1 1 that is compensated for entrained gas/air in the liquid or 
mixture 12. The apparatus 300 includes a device 302 for measuring the gas volume 
fraction of the mixture, which is substantially similar to that described herein before, and 

20 an electromagnetic flow meter 304 (also know as a Magmeter) for measuring the 

uncompensated volumetric flow rate of the mixture 12 within the pipe 14. As shown, the 
GVF device 302 provides a signal 8 indicative of the gas volume fraction of the mixture 12 
to the electromagnetic flow meter 304. The magmeter 304 then compensates (or corrects) 
volumetric flow measurement 308 for entrained gas/air in response to the gas volume 

25 fraction signal 8. For example, this can be accomplished in accordance with the following 
equation: 

Compensated Vol. Flow Rate = Vol. Flow Rate (1 - Gas Volume Fraction) 
The electromagnetic flow meter 308 may comprise an electromagnetic flow meter as 
described hereinbefore, such as that 8700 Series Magmeter manufactured by Rosemount. 
30 However, one will appreciate that the entrained air meter portion may used to compensate 
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or correct any volumetric flow meter that is able to provide a volumetric flow meter 
measurement. 

While the gas volume fraction signal 8 was provided to the magmeter 304 to 
compensate the volumetric flow rate, the invention contemplates that the correction for 
5 entrained gas/air may be performed in the GVF device 302. 

Similarly, this GVF device 302 for measuring entrained air or gas volume fraction 
measurement GVFair 8 may be used with any consistency meter to correct for errors 
introduced into a measurement by entrained gas/air. For instance, a consistency meter will 
show an error when entrained air exists in the mixture. The present invention may be used 
1 0 to correct for this error. 

Referring to 14, the present invention contemplates an apparatus 400 for providing a 
consistency measurement 1 1 that is compensated for entrained gas/air in the mixture or 
slurry 12 (e.g. paper/pulp slurry). The apparatus 300 includes a device 302 for measuring 
the gas volume fraction of the slurry 12, which is substantially similar to that described 
1 5 herein before, and a consistency meter 402 for measuring the uncompensated consistency 
measurement of the mixture 12 within the pipe 14. As shown, the GVF device 302 
provides a signal 8 indicative of the gas volume fraction of the slurry 12 to the consistency 
meter 402. The consistency 402 then compensates (or corrects) consistency measurement 
404 for entrained gas/air in response to the gas volume fraction signal 8. For example, this 
20 can be accomplished in accordance with the following equation: 

Cmixture = Uncompensated Consistency Measurement - (1 .4)(GVFair) 
The consistency meter 402 may comprise a microwave consistency meter as 
described hereinbefore, such as that 8700 Series Magmeter manufactured by Rosemount. 
However, one will appreciate that the consistency meter portion may used to compensate or 
25 correct any consistency meter that is able to provide a consistency measurement. 

While the gas volume fraction signal 8 was provided to the consistency meter 402 to 
compensate the consistency measurement as shown in Figs. 14, the invention contemplates 
that the correction for entrained gas/air may be performed in the GVF device 302. 
Referring to Fig. 15, the present invention contemplates the consistency meter 402 in 
30 accordance with the present invention receives and processes the measured pressure input 
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6 and the gas volume fraction input 8 (compensated for pressure) to provide a corrected 
consistency output 406 to the plant DCS 410. 

In Fig. 16, the plant DCS 410 in accordance with the present invention receives and 
process the uncorrected consistency measurement 404 from the consistency meter 402 and 
5 a GVF input 8 (compensated for pressure) and provides a corrected consistency output 
406. In Fig. 17, the plant DCS 410 in accordance with the present invention receives and 
processes an uncorrected consistency measurement, a pressure measurement, and GVF 
measurement and provides a corrected consistency measurement 406 and a GVF 
measurement (compensated for process pressure). 

10 Figs. 1 8 - 20 show data from an apparatus 400 in accordance with the present 

invention that measured and corrected for entrained air within a slurry propagating through 
a pipe. Fig. 18 shows the measured consistency from a consistency meter 402, the 
measured GVFair for a GVF meter 302 and the corrected consistency 406 over an eight 
hour time period. Fig. 19 shows the measured consistency from a consistency meter 402, 

1 5 the measured GVFair for a GVF meter 302 and the corrected consistency 406 for a slurry 
having macro-bubbles flowing therein. Fig. 20 shows the measured consistency from a 
consistency meter 402, the measured GVFair for a GVF meter 302 and the corrected 
consistency 406 for a slurry having micro-bubbles flowing therein. Consequently, as 
shown, the present invention eliminates the big bubble/small bubble analysis/compensation 

20 that other flow meters attempt to do to compensation for entrained gas/air in the process 
flow 12. 

It should be understood that any of the features, characteristics, alternatives or 
modifications described regarding a particular embodiment herein may also be applied, 
used, or incorporated with any other embodiment described herein. 
25 Although the invention has been described and illustrated with respect to exemplary 

embodiments thereof, the foregoing and various other additions and omissions may be made 
therein and thereto without departing from the spirit and scope of the present invention. 
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